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Abstract. We investigate the fine structure of the sunspot penumbra by means of a model that allows for a flux 
tube in horizontal pressure balance with the magnetic background atmosphere in which it is embedded. We apply 
this model to spectropolarimetric observations of two neutral iron lines at 1.56 fim and invert several radial cuts 
in the penumbra of the same sunspot at two different heliocentric angles. In the inner part of the penumbra we 
find hot flux tubes that are somewhat inclined to the horizontal. They become gradually more horizontal and 
cooler with increasing radial distance. This is accompanied by an increase in the velocity of the plasma and a 
decrease of the gas pressure difference between flux tube and the background component. At large radial distances 
the flow speed exceeds the critical speed and evidence is found for the formation of a shock front. These results 
are in good agreement with simulations of the penumbral fine structure and provide strong support for the siphon 
flow as the physical mechanism driving the Evershed flow. 
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1. Introduction 

The picture of the magnetic fine structure of the penumbra 
has strongly evolved over the last decade (e.g., Degenhardt 
& Wiehr 1991; Title et al. 1993; Solanki & Montavon 1993; 
Rinmrele 1995; Stanchfield et al. 1997; Westendorp Plaza 
et al. 1997; Schlichenmaier et al. 1998a; Scharmer et al. 
2002; see Solanki 2003, Bellot Rubio 2003 and Thomas & 
Weiss 2004 for an overview). It is now accepted that the 
penumbral magnetic field is uncombed , i.e. inclined at least 
in two different directions on a small scale. There is also 
considerable evidence that the more horizontal component 
must be in the form of flux tubes, although the diameter 
of these flux tubes is still a matter of debate (Sanchez 
Almeida 1998, 2001; Martinez Pillet 2000, 2001). These 
flux tubes carry the Evershed flow (Evershed 1909; Title 
et al. 1993; Westendorp Plaza et al. 2001a,2001b; Bellot 
Rubio et al. 2003,2004; Borrero et al. 2004). Many of the 
tubes return to the solar interior within the penumbra 
(Westendorp Plaza et al. 1997; del Toro Iniesta et al. 2001; 
Mathew et al. 2003; Borrero et al. 2004) and along with the 
magnetic flux a large fraction of the mass flux carried by 
the Evershed flow returns also to the solar interior within 
the penumbra or just outside it (Borner & Kneer 1992; 
Solanki et al. 1994, 1999). 


* Present address: High Altitude Observatory, 3450 Mitchell 
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The combination of these results raises questions re¬ 
garding the commonly considered physical mechanism 
driving the Evershed flow. Since the wave hypothesis is 
ruled out (Biinte & Solanki 1995) and episodic Evershed 
flow produced when a flux tube falls and drains (Wentzel 
1992) faces difficulties due to the relative immutability of 
the penumbral magnetic pattern (Solanki & Riiedi 2003), 
the most widely accepted mechanism is that the flow is 
caused by a gas pressure gradient between the upflow- 
ing and the downflowing footpoints, also referred to as 
the siphon flow mechanism (Meyer & Schmidt 1968; see 
also Thomas 1988). Commonly, this pressure gradient is 
thought to be produced by a difference in the field strength 
between the footpoints which, due to horizontal pressure 
balance, leads to a gas pressure difference (e.g. Degenhardt 
1989,1991; Thomas & Montesinos 1991,1993, Montesinos 
& Thomas 1997). All else being equal, the gas flows from 
the footpoint with lower field strength to the footpoint 
with the higher field strength. However, if most of the gas 
flows only within the penumbra, then due to the roughly 
factor of 2 larger magnetic field strength at the inner 
boundary of the penumbra compared to the outer edge, 
one would naively expect the gas to flow inwards, contra¬ 
dicting observations. Montesinos & Thomas (1997) have 
argued that this radial decrease of the field strength is 
only apparent, being caused by different r = 1 levels at 
the footpoints and the fact that local intense magnetic flux 
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concentrations at the outer penumbral edge could not be 
easily resolved. 

A possible resolution of this dilemma was noticed by 
Borrero et al (2004; hereafter Paper I), who found that 
whereas the strength of the inclined magnetic compo¬ 
nent drops very rapidly in the radial direction, as re¬ 
quired by a global magnetic structure close to potential 
(see Jahn 1989), the horizontal component carrying the 
Evershed flow shows far less variation (cf. Riiedi et al. 
1998, 1999). The analysis in Paper I was incomplete in the 
sense that the two components were independent of each 
other. Here we overcome this shortcoming and take into 
account the flux-tube structure of the field and the pres¬ 
sure balance between components. We apply a powerful 
inversion technique (described in Section 2) to spectropo- 
larimetric observations of infrared Zeeman sensitive spec¬ 
tral lines (Section 3). General results are presented and 
discussed in Section 4. In Section 5 we discuss the genera¬ 
tion of Net Circular Polarization by the uncombed model 
and the implications for the typical size of the penumbral 
flux tubes. Section 6 describes our results in the frame¬ 
work of the theoretical models employed to explain the 
Evershed effect in terms of gas pressure differences. Our 
main findings and conclusions are summarized in Section 
7. 


2. Description of the penumbral model and 

Stokes profile inversion 

The analysis of spectropolarimetric observations of the 
sunspot penumbra by means of Stokes profile inversion 
has, so far, either considered 2 distinct components (as¬ 
suming the physical quantities to be constant with depth: 
Bellot Rubio et al. 2003,2004) or 1 component (allowing 
gradients to be present; Westendorp Plaza et al. 1997, 
2001a, 2001b; Bellot Rubio et al. 2002; Mathew et al. 
2003). In Paper I we suggested that the uncombed penum¬ 
bral model described by Solanki & Montavon (1993) pro¬ 
vides a picture for the penumbral fine structure that is 
able to encompass the results of these investigations. This 
conclusion was based on the application of a Stokes profile 
inversion technique assuming the two different geometries 
mentioned above, together with considerations on how the 
area asymmetry of the circular polarization profiles is in¬ 
fluenced by gradients in the magnetic and kinematic strat¬ 
ifications. 

In this paper we carry out Stokes profile inversions 
based on the uncombed model. This consists of a flux 
tube embedded in a magnetic surrounding atmosphere. 
The forward modelling for the considered geometry has al¬ 
ready been addressed (see e.g. Degenhardt & Kneer 1992, 
Solanki & Montavon 1993; Martinez Pillet 2000). The ba¬ 
sic geometry is illustrated in Figure 1. The simplest repre¬ 
sentation of the uncombed model is in terms of two rays: 
the first ray passes along the surrounding atmosphere (ver¬ 
tical dot-dashed line pointing towards the observer) while 
the other ray cuts both surrounding atmosphere and the 


OBSERVER 



Fig. 1. Example of the geometrical scheme used in the 
inversion. The radiative transfer equation is solved along 
the 2 rays (dashed and clot-dashed lines) representing the 
embedded flux tube and surrounding atmosphere respec¬ 
tively. 7 S and refer to the inclination of the magnetic 
field vector with respect to the observer. In this picture 
for simplicity the heliocentric angle is 6 = 0 and q t = 90°. 

flux tube (vertical dashed line) 1 . Let us denote with % s 
any of the magnetic and kinematic physical quantities for 
the surrounding atmosphere 2 . For the flux tube compo¬ 
nent we adopt the following form for x t (z): 

v _ f Xt if z€ [z 0 —R t ,z 0 + R. t ] m 

^ \ Xs otherwise 1 

where Xt on the right hand side of Eq. 1 is height in¬ 
dependent, Zq is the height where the axis of the flux 
tube is located and R t is its radius. Note that Eq. 1 im¬ 
plies that physical stratifications along the flux tube at¬ 
mosphere (vertical dashed line in Fig. 1) are the same as 
in the surrounding atmosphere above and beneath the flux 
tube: z < z Q — Rt and z > z 0 + Rt . At the flux tube lower 
and upper boundaries (z* = zq ± Rt) the physical quan¬ 
tities suffer a jump whose magnitude is = xt — Xs- 
These jumps/gradients are the essential ingredients to ex¬ 
plain the Net Circular Polarization (NCP) observed in 
the sunspot penumbra (Solanki & Montavon 1993; cf. 
Martinez Pillet 2000). 

The inversions based on this geometry have been car¬ 
ried out using the inversion code SPINOR (see Frutiger 
et al. 1999; Frutiger 2000). The code performs spec¬ 
tral line synthesis in Local Thermodynamic Equilibrium 
(LTE) and employs a Levenberg-Marquart nonlinear least- 
squares x 2 minimization algorithm (Press et al. 1986), 
whereby derivatives are calculated through numerical re¬ 
sponse functions (RFs; see Ruiz Cobo & del Toro Iniesta 
1992). The free parameters allowed in the inversion are, 
for the surrounding atmosphere: Vlos.s (line of sight ve- 

1 Note that this geometrical simplification in terms of two 
rays implies that we are considering the cross section of the 
flux tube to be square instead of circular as in Fig. 1. 

2 Xs is assumed to be height independent. 
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locity), B s (magnetic field strength), (magnetic field in¬ 
clination with respect to the observer), <p s (magnetic field 
azimuth) and T s (t5 = 1) for temperature, where 75 is the 
optical depth at a reference wavelength of 5000 A. The 
height dependence of T s (r) is taken from the penumbral 
model by del Toro Iniesta et al. (1994). When T s (ts = 1) is 
changed by AT the temperature at all heights is changed 
additively: T s (r) = T Sj 0 ld (r) + AT. 

For the flux tube the free parameters are: VLOS.t, Bt , 
7 t, <f>t, macro and microturbulent velocities, u maC) t and 
i> m ic,t) and T t (r 5 = 1) (where again T t (r) is taken from 
the mean penumbra model of del Toro Iniesta et al.). In 
addition, zo, Rt and at (fractional horizontal area cov¬ 
ered by the flux tube component with respect to the total 
area) are also allowed to change. Finally we employ a stray 
light correction to model the contribution of light from 
the neighbouring granulation that enters into the spectro¬ 
graph’s slit. To this end we used the two component model 
for the quiet sun from Borrero & Bellot Rubio (2002) to 
produce synthetic intensity profiles, I q , of the observed 
spectral lines (see Sect. 3) and combined it with the emer¬ 
gent spectrum of the pure penumbral profiles using a filling 
factor a q which is also retrieved from the inversion. This 
results in a total of 16 free parameters. This compares pos¬ 
itively with the inversions carried out in Paper I, which 
used a total number of 23 and 18 free parameters for the 
1 component and 2 component inversions, respectively. 

The radiative transfer equation is integrated using the 
Hermitian Approach (see Bellot Rubio et al. 1998) for 
each ray separately. The Stokes profiles from quiet sun, 
flux tube and surrounding atmosphere are finally com¬ 
bined using the filling factors at and a q : 


S( A) = a q S q ( A) + (1 - a q )[(l - a t )S s ( A) + a t 5 t (A)] ( 2 ,) 

where S represents the Stokes vector (I,Q,U,V): as aris¬ 
ing from the surrounding atmosphere (dot-dashed ray in 
Fig. 1) S s , from the ray cutting the flux tube (dashed 
ray in Fig. 1) St, as well as the quiet sun contribution 
S q = (lq, 0 , 0 , 0 ) which is a non polarized contribution and 
only affects the total intensity profiles. These synthetic 
profiles are now compared with the observations and the 
free parameters are changed according to the RF until the 
minimum x 2 is reached. A detailed study of the numerical 
performance of this procedure as well as the uniqueness of 
the retrieved atmosphere under different levels of noise in 
the observations is presented by Borrero et al. (2003a). 

The radiative transfer is always performed in the op¬ 
tical depth scale, but the flux tube and its force balance 
with its surroundings are more naturally described in the 
geometrical depth scale. The correct functioning of this 
interplay requires a sufficiently complex procedure. First, 
a geometrical height scale is assigned to the surrounding 
atmosphere following the strategy outlined in Gray (1992) 


and integrating the hydrostatic equilibrium equation (as¬ 
suming a force free situation V x B |j B), 

^ = ( 3) 

OT s K c>s 

where g is the solar surface gravitational acceleration 
and r is the optical depth computed using the contin¬ 
uum opacity, k c , at 5000 A. This requires an estimate 
of the gas pressure at the top of the tabulated atmo¬ 
sphere: Pg tS {T s ^ max ). With this, the gas pressure stratifica¬ 
tion Pg tS {T s ) is obtained in the surrounding atmosphere, 
and since the temperature is obtained from the inversion, 
the equation of state (ideal gas law including a variable 
mean molecular weight to account for the ionization of 
the different species) provides the density: p s (r s ). The re¬ 
lation dr s = —p s K c , s dz is now integrated setting z = 0 at 
log 75 = 0 and thus defining the geometrical height scale. 

The gas pressure in the tube component is obtained 
under the assumption of total pressure balance with the 
surroundings at the height of the axis of the flux tube, 

P t ( z ) = / ^s,s( z ) + ) if 2 G [zq - R t , z 0 + -R t ] M-v 

s ’ y Pg,s{z) otherwise 

This assumption is valid as long as the magnetic field 
of the external atmosphere does not penetrate into the 
flux tube and vice versa (Kippenhahn & Mollenhof 1975, 
Chap.3). Since P g j{z) is now known and Tt(z) was ob¬ 
tained from the inversion, the density pt{z) can be evalu¬ 
ated by using the idal gas equation, and thus a new optical 
depth scale for the atmosphere containing the flux tube 
can be obtained through the relation: dr t = — ptn c ,tdz. 
For the integration of this last equation a boundary con¬ 
dition is employed which implies that for z > zo + Rt 
the surrounding and flux tube components must have the 
same z values: z(r t ) = z(t s ). Note that the obtained flux 
tube density, pt(z), does not satisfy vertical hydrostatic 
equilibrium. 

3. Observations and data reduction 

The active region NOAA AR 8706 was observed on 
Sep 21st, 1999 and Sep 27th, 1999 at /./ = cosd = 
0.51 and p = 0.91 respectively, using TIP (Tenerife 
Infrared Polarimeter, Martinez Pillet et al. 1999) attached 
to the spectrograph of the 70cm German VTT of the 
Observatorio del Teide. Here, 9 is the heliocentric angle. 
The recorded spectral region contains the full polarization 
profiles of the pair of Fe I lines 15648.5 A (g = 3) and Fe I 
15652.8 A(g e ff = 1.53). The wavelength sampling is about 
29 mA. The diagnostic properties of these lines have been 
discussed by Solanki et al. (1992). They are formed in 
the deep photosphere as a result of their high excitation 
potentials and the low continuum opacity at this wave¬ 
length. These lines sample a relatively narrow layer not 
wider than log 75 = [0.5, —2] (see Bellot Rubio et al. 2000; 
Mathew et al. 2003). In the umbra the second neutral 
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iron line is heavily blended by two molecular OH lines at 
15651.9 and 15653.5 A. These lines are Zeeman sensitive 
and belong to the Meinel system 3,1 (see Berdyugina & 
Solanki 2002 for details). Their equivalent widths greatly 
decrease towards the penumbra. However they can still be 
seen clearly at the umbra-penumbra boundary. In order 
to accurately analyze both atomic lines we also compute 
the OH lines. The amployed atomic parameters for the ob¬ 
served Fe I lines were taken from Borrero et al. (2003b), 
while for the OH lines the values given by Abrams et al. 
(1994) and Mies (1974) were used (see Table 1). 

The usual data reduction procedures for TIP were in 
general followed; we proceeded carefully at several points, 
however. Firstly the neutral iron line, Fe I 15648.5, ap¬ 
pears to be blended by a telluric line in its red wing. This 
complication was overcome by fitting the average quiet 
Sun profile of the second line with an appropriate model 
atmosphere. This model was then used to synthetize the 
first iron line. The ratio between the average observed and 
computed line allowed us to recover the shape of the tel¬ 
luric blend, which was subsequently eliminated from the 
rest of the profiles. This procedure, although not perfect, 
does not introduce modifications to the original profile 
(equivalent width and line core intensity) larger than ~3 
%. 

Secondly, the continuum correction posed a consider¬ 
able problem. The four quadrant configuration of the TIP 
camera produce small gradients in the continuum inten¬ 
sity that remain after applying the flat field correction. 
To account for this we compared our spatially averaged 
flat field with the infrared FTS atlas from Livingston & 
Wallace (1991) and defined several wavelength positions 
where the continuum should be reached. A smooth sec¬ 
ond order polynomial is interpolated over those selected 
points and used to bring the continuum in the flat field 
image to the same level as the FTS continuum. Thirdly, 
the polarization signals were corrected for residual cross 
talk using the statistical approach described in Collados 
(2001; see also Schlichenmaier & Collados 2002). Finally, 
the wavelength was calibrated by assigning the laboratory 
wavelengths to the cores of the Fe I lines in the average 
quiet Sun profiles. To account for the effects of the gran¬ 
ulation we shifted this wavelength scale by —400 nr s -1 , 
which is the approximate value for the convective blueshift 
as deduced from the quiet sun model of Borrero & Bellot 
Rubio (2002). Second order corrections due to the different 
viewing angles (Balthasar 1985) have not been performed. 

The seeing conditions were rather good during the ob¬ 
servations, with the granulation being clearly discernible 
in the reconstructed continuum images (see Figure 2; left 
panels). By calculating the power spectrum of the contin¬ 
uum intensity in the neighbouring quiet Sun we estimate 
the spatial resolution to be about 1 arc sec. In Fig. 2 (right 
panels) the maps of the total circular polarization are also 
shown. Due to projection effects the polarity inversion line 
(region in the limb side of the penumbra where the Stokes 
V signal changes its sign, i.e., where the average magnetic 
field is perpendicular to the observer) appears in the outer 


Table 1 . Atomic and molecular parameters of the ob¬ 
served spectral lines. Ao represents the laboratory central 
wavelength, \\ the excitation potential of the lower energy 
level, and log gf the logarithm of the oscillator strength 
times the multiplicity of the level. The parameters a and 
a (in units of Bohr’s radius, ao) are used to calculate the 
broadening of the lines by collisions with neutral hydro¬ 
gen atoms as resulting from the ABO theory (Barklem & 
O’Mara 1997). The last column gives the effective Lande 
factor of the transition, g e ff. For the molecular lines Itj, Jl, 
Vjj and Vl represent the upper/low multiplets sublevels 
and vibrational levels respectively. Jl stands for the ro¬ 
tational number of the lower level. Finally, the oscillator 
strength of the molecular transition is given. 


Atom Ao Xi l°g gf a cr g eS 

_ (A ) (eV) (dex) _(ao)_ 

FeT 15648.515 5.426 -0.675 0.229 977 TOO 

Fe I 15652.874 6.246 -0.043 0.330 1444 1.53 

Molecule Ao Branch Vu-Vl Jl / 

OH 15651.895 P Tl 3A T5 0.8xl(T o 

OH 15653.478 P 1-1 3-1 6.5 O.SxlO -6 


penumbra when the sunspot is observed near disk center 
(n = 0.91) but lies closer to the umbra when the spot is 
located closer to the limb (p = 0.51). Note that the arrows 
in Fig. 2 point towards the solar disk center. Part of this 
data set has been analyzed in Paper I (see also Mathew 
et al. 2003, 2004). 

4. Results and discussion 

We have inverted individually all the pixels along the ra¬ 
dial cuts shown in Fig. 2 (5 for each heliocentric angle). 
Each cut contains approximatively 20 pixels and ranges 
roughly from r/R = 0.4 to r/R = 0.9 (where R is the 
penumbral radius, indicated by the external contour in 
Fig. 2). We have chosen the cuts such that they all lie 
on the limbward side of the penumbra and near the line 
of symmetry (i.e. the line connecting the sunspot’s center 
and the center of the solar disk), since this is where the 
flux tubes leave the most distinctive fingerprints on the 
observed profiles, thus allowing for a reliable determina¬ 
tion of their properties, as shown in Paper I. 

4.1. Example 

In Fig. 3 we present an example of the observed (filled 
circles) and fitted (solid lines) circular polarization pro¬ 
files for a penumbral point. These multi-lobed profiles can 
be successfully fitted by adding together two Stokes V 
profiles. The profile from the surrounding atmopshere, 14, 
is produced by a positive polarity magnetic field (7 < 
90°; dot-dashed line) with zero line-of-sight velocity (zero 
crossing is at Ao), and the Stokes V profile from the ray 
cutting through the flux tube, 14 , which is produced by a 



Borrero et al: On the nature of the Evershed flow 


5 



arc sec 

Continuum Intensity at 1.56 /irr 



10 20 30 40 


arc sec 


arc sec 

Total Circular Polarization 



arc sec 


Fig. 2. Top panels: NOAA 8706 observed in 21 September 1999 at a heliocentric angle p = 0.51 (left: continuum 
intensity map at 1.56 pm; right: total circular polarization map for Fe I 15648.5 A). Bottom panels: NOAA 8706 
observed in 27 September 1999 at a heliocentric angle /i=0.91 (left: continuum intensity map at 1.56 /im; right: total 
circular polarization map for Fe I 15648.5 A). The arrows point towards the direction of the solar disk center. The two 
inner-most contours enclose the umbral-penumbral boundary. The external contour defines the sunspot radius r = R, 
at each position angle. These three contours have been defined as 0.45/ c ,0.65/ c and 0.85J c , where I c represents the 
average continuum intensity of the quiet Sun. The radial cuts selected for our analysis are also shown. 


negative polarity magnetic field (7 > 90°; dashed line) car¬ 
rying a flow directed away from the observer (zero crossing 
is red shifted with respect to Ao). The filling factor of the 
flux tube atmosphere a t and of the stray light contribution 
a q are applied to V s and Vt to obtain the final emergent 
profile (solid line) according to Eq. 2. 


In Fig. 4 we present the atmospheric stratifications re¬ 
sulting from the inversion of the profiles shown in Fig. 3 
using the uncombed model (see Sect. 2). Outside the flux 
tube the magnetic field strength is about 2000 Gauss and 
is inclined with respect to the observer by about 75°. 
Inside the flux tube the magnetic field is weaker (~ 1300 
Gauss) and more inclined (7 ~ 100°). While the surround¬ 
ing atmosphere is basically at rest, Vlos.s — 0.2 km s _1 , 
we detect red shifts of IAos,t — 1.1 km s _1 inside the 
flux tube. In addition, the results from the inversion indi¬ 
cate that the flux tube is hotter, by roughly 500 K, than 
the surrounding atmosphere. The presence of the discon¬ 
tinuity along the flux tube atmosphere produces the area 
asymmetry in Stokes V, AAfit, marked in Fig. 3. 


4.2. General properties 

We have taken the individual results of the inversions of 
the considered radial cuts (separately for each heliocentric 
angle) at the geometrical height z = Zq (i.e. at the location 
of the tube axis), and calculated the averages at each ra¬ 
dial position in the sunspot. This is done individually for 
the flux tube component and the surrounding magnetic 
field. The results are presented in Fig. 5 and 6 for the 
temperature, line-of-sight velocity, the flux tube’s filling 
factor at, magnetic field zenith angle and magnetic field 
strength. In order to compare the magnetic field inclina¬ 
tion deduced for sunspots located at different heliocentric 
positions we have converted from the observer’s reference 
frame, 7 , to the local reference frame. Therefore, we plot 
now the zenith angle, £ (Fig. 6 ; top panels). A zenith an¬ 
gle smaller than, equal to or larger than 90° indicates that 
the magnetic field is inclined upwards, lies parallel to or 
is inclined downwards with respect to the solar surface. 
All in all, the results for the two sets of observations are 
remarkably similar in spite of the different viewing angles 
and the time difference of 6 days between the two obser¬ 
vations. 
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Fig. 3. Example of the observed (dots) and fitted (solid lines) Stokes V profiles (left: 15648.5 A; right: 15652.5 A) 
for a penumbral point. The fitted profile is obtained by linearly combining the profile emerging from the surrounding 
component (dot-dashed line) and the profile emerging from the ray piercing the flux tube (dashed line). The employed 
filling factors, at and a q , as well as the area asymmetry of the observed and fitted profiles are also given. 
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Fig. 4. Temperature (top-left panel), magnetic field strength (top-right panel), magnetic field inclination (bottom-left) 
and line-of-sight velocity (bottom-right) for the flux tube atmosphere (dashed lines) and its surroundings (dot-dashed 
lines) as a function of the geometrical depth, obtained from the inversion of the profiles in Fig. 3. The flux tube’s 
radius returned by the inversion is 125 km and its central position is zq = 150 km. 


In the inner penumbra, we detect nearly, but not com¬ 
pletely horizontal flux tubes (£ t — 70—80°) that are hotter 
than their surroundings by about 500-1000 K. These flux 
tubes carry most of the Evershed flow, with LOS velocities 
in the inner penumbra ranging from 0.5 to 3 km s -1 . The 
magnetic field strength in the flux tubes is around 1500 


Gauss. The atmosphere surrounding these flux tubes pos¬ 
sesses a more vertical (£ s ~ 20 — 40°) and stronger mag¬ 
netic field (B ~ 2300-2500). No signatures of the Evershed 
flow are detected here. 

As the radial distance increases, the flux tubes cool 
down to temperatures similar to those of the surroundings. 
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At large radial distances r/R > 0.7 the temperature in the 
flux tube component decreases even below the surrounding 
temperature, although only slightly (~ 200-300 K). At the 
same time the tubes become more horizontal, reaching 
£t — 90°, and point slightly downwards with respect to 
the solar surface, — 95 — 100°, near the outer edge of the 
penumbra, i.e. at r/R > 0.8. In addition, their magnetic 
field strength decreases slowly to 800-1000 Gauss at r/R = 
0.9 while the LOS velocity increases monotonically in both 
spots, although for fi = 0.51 it suffers a sudden drop near 
the outer penumbral boundary. The filling factor of the 
flux tubes (Fig. 5; bottom panels) increases continuously 
from very small values at the inner penumbra a t — 0.2 — 
0.3 until they cover almost all the resolution element at the 
outer boundary a t — 0.9. This can be interpreted either 
as an increase in the horizontal cross section of the flux 
tubes or as an increase in the number of flux tubes per 
resolution element. 

The surrounding atmosphere exhibits a quite differ¬ 
ent behaviour. The magnetic field strength decreases very 
rapidly towards the outer penumbra, reaching similar val¬ 
ues to those of the flux tube’s magnetic field strength (~ 
800-1000 Gauss) at r/R > 0.8. The inclination of the 
magnetic field increases slightly with radial distance up to 
</ s ~ 40 — 50°. The LOS velocities remain small throughout 
the penumbra. 

Before discussing the individual results in more detail 
there are several points which need to be clarified. First 
of all, we want to stress that the results obtained from the 
inversions of the radial cuts when the spot is near the disk 
center (Figs. 5 and 6; left panels) are more reliable at in¬ 
termediate to large radial distances, while results inferred 
from the sunspot at large heliocentric angles (Figs. 5 and 
6 ; right panels) are more reliable in the inner penumbra. 
The reason for this is that the magnetic neutral line is 
located in the outer penumbral for sunspots near the disk 
center, /r = 0.91, but it shifts towards the umbra as the 
sunspot is located closer to the limb (/j = 0.51). The mag¬ 
netic neutral line is where multi-lobed Stokes V profiles 
are commonly observed (Sanchez Almeida & Lites 1992; 
Sclilichenmaier & Collados 2002). As already discussed in 
Paper I inferred parameters from the inversion are more 
reliable for exactly such complex profiles, since the sig¬ 
nature of the two unresolved components (flux tubes and 
magnetic surrounding) can be better distinguished there 
(e.g. Fig. 3), allowing for a more reliable inference of their 
properties. In general, however, larger and more accurate 
line-of-sight velocities are obtained from the sunspot near 
the limb (at all radial distances) since the Evershed flow 
is more aligned with the line-of-sight there. 

We note that the possible size of the flux tubes can 
be much smaller than our spatial resolution of 1 arc sec 
(Sutterlin 2001; Sutterlin et al. 2004; Scharmer et al. 
2002; Rouppe van der Voort et al. 2004). Threfore the de¬ 
duced properties are unlikely to correspond to one single 
flux tube but rather to some average over all the possi¬ 
ble penumbral fibrils contained in the resolution element. 
However, as one can see from Fig. 5 and 6, the inferred 


properties are similar to those expected for a single flux 
tube which crosses the penumbra from its inner to the 
outer boundary. 


5. Width of the penumbral filaments and Stokes 
V area asymmetry 

The area asymmetry, 8A, of the circular polarization is 
defined as: 


fV( A)dA 
/|V(A)|dA 


( 5 ) 


It is different from zero whenever gradients along the line- 
of-sight of the magnetic field vector and LOS velocity are 
present (Landolfi & Landi degl’Innocenti 1996). Solanki & 
Montavon (1993) realized that the huge gradients needed 
to reproduce the area asymmetry observed in the sunspot 
penumbra with the visible Fe I lines at 6301.5 and 6302.5 
A (Sanchez Almeida & Lites 1992) could be interpreted 
as a horizonal flux tube embedded in a more vertical back¬ 
ground. As the line of sight crosses the tube’s boundaries 
the physical stratifications describing the atmosphere suf¬ 
fer a jump that is directly responsible for the generation 
of the area asymmetry in Stokes V. Schlichenmaier at al. 
(2002) and Muller et al. (2002) further investigated this 
issue and pointed out that the area asymmetry observed 
in the visible Fe I 6301 A lines is dominated by jumps in 
the magnetic field inclination, while the area asymmetry 
observed in the infrared 1.56 /cm lines can be explained in 
terms of jumps in the magnetic field azimuth. Both cases 
are compatible with the uncombed penumbra proposed by 
Solanki & Montavon (1993). 

As already described in Sect. 2, in the uncombed 
model, the discontinuities needed to produce asymmet¬ 
ric circular polarization profiles are located at the tube’s 
boundaries: z* = Zo±Rt, with R t and Zq (the tube’s radius 
and central position) being free parameters of the inver¬ 
sion. The radius we obtain from the inversion of all the 
pixels contained in the 10 radial cuts in Fig. 1 in the same 
sunspot at two different heliocentric angles is, consistently, 
125 Kilometres, which is precisely the maximum value we 
allow for R t 3 . This means that the inversion tries to make 
the flux tube as thick as possible. The area asymmetry of 
the synthetic profiles, Mfit only poorly reproduces the 
observed one: Mobs ( see Fig. 7; top panels). In most cases 
the uncombed model produces an area asymmetry which 
is very little or almost zero, and varies over a smaller range 
than Mobs- 

The small observed area asymmetries (average of 
Mobs ~ 3 %) clearly indicate that SA is a parameter 
which plays a minor role in the inversion since the shape 
of the circular polarization profiles can be successfully fit¬ 
ted by means of profiles showing little or even no area 

3 This maximum value is set in the inversion code so that 
the upper boundary of the flux tube lies below the top of the 
tabulated atmosphere. This ensures that the condition z(rt) = 
z{t s ) for z > zo T Rt (see Sect. 2) can be applied. 
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Fig. 5. Top panels: temperature difference between the flux tube atmosphere and its surroundings as a function of 
radial distance from the spot center. Middle panels: radial variation of the line-of-sight velocities inside (dashed lines) 
and outside the flux tubes (dashed-dotted lines). Bottom panels: radial variation of the flux tube’s filling factor, at- 
Left panels: sunspot at // = 0.91. Right Panels: sunspot at /i = 0.51. All quantities refer to the central position of the 
tube, z = zq. Shaded areas denote the maximum deviations from the average at each radial position, obtained from 
the individual inversions of all the radial cuts. The arrow indicates the approximate radial position where the flux 
tube experiences a final temperature enhancement (see Sect. 6.4). 


asymmetry at all (as the 2C model in Paper I). Even ex¬ 
tremely strange profiles, such as those presented in Fig. 3, 
can be fitted reasonably well. This is because the circular 
polarization profiles of the spectral lines used in this work 
are mainly affected by the presence of two different polar¬ 
ities in the resolution element (background with 7 < 90° 
and flux tube with 7 > 90°) rather than by discontinuities 
along the line of sight in the physical parameters. This is 
due to the large Zeeman splitting of these loines and the 


fact that they are not strongly saturated (see Grossmann- 
Doerth et al. 1989). Note that Borrero et al. (2004; paper 
I) also reached the same conclusion and suggested that 
using the visible iron lines at Fe I 6301 A, whose area 
asymmetry is far more sensitive to such discontinuities, 
might help to further constrain the size of the penumbral 
filaments. For the case of the infrared Fe I lines at 1.56 
/jm they constrained a small region where the area asym- 
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Fig. 6. The same as Fig. 5, but now for the radial variation of the magnetic field zenith angle (top panels) and magnetic 
field strength (bottom panels). 


metry of these lines is sensitive to gradients along the line 
of sight: logr 5 € [0, —0.5]. 

In Fig. 7 (middle panels) the position (in the opti¬ 
cal depth scale) of the tube’s lower boundary is plotted 
as a function of r/R. At r/R < 0.7 the lower boundary 
is located at logrs > 0. At the outer penumbra, how¬ 
ever, the lower boundary shifts to higher layers: logrs ~ 
—0.25, where the discontinuity is effective in generating 
area asymmetry (see Paper I). Note that for most of the 
penumbra the lower boundary also remains below the 
logn.se level (dashed line). For comparison we plot in 
Fig. 7 (bottom panel) the observed, <L4 obs (solid line), 
and fitted (dashed line), 5 Afit, area asymmetry for all 
the considered radial cuts versus r/R. <L4obs increases ra¬ 
dially in the penumbra, from small negative values (~ —2 
%) up to larger positive ones: ~ 7 % . Interestingly the 
fitted area asymmetry displays a similar behaviour, spe¬ 
cially in the outer penumbra, where observed and fitted 
values become similar. 

The combination of these results leads us to conclude 
that, in the inner-intermediate penumbra, where the ob¬ 
served area asymmetry in the circular polarization profiles 
of the Fe I 1.56 pm lines is small, discontinuities are not 
important to reproduce the profiles and therefore the tube 
radius is a parameter which is not well constrained from 
the inversion. However, in the outer penumbra, the ob¬ 
served area asymmetry becomes large enough to turn into 


an important ingredient to successfully fit the circular po¬ 
larization profiles. In these regions, the lower boundary of 
the flux tube is located at heights where it is effective in 
generating area asymmetry. As the upper boundary re¬ 
mains unseen we still cannot draw any reliable conclusion 
on the actual size of the penumbral filaments. 

Given that Fe I lines at 1.56 pm do not see much of 
the tube’s boundaries, our uncombed model could in prin¬ 
ciple be simplified into a two component model where all 
physical quantities are constant with height (e.g., mak¬ 
ing the tube’s radius infinite in Fig. 4). Such models have 
been previously used to study the fine structure of the 
penumbra (Bellot Rubio et al. 2003,2004; Borrero et al. 
2004). However a feature included in the uncombed model 
but has been neglected by previous two component mod¬ 
els remains: the use of the total pressure balance between 
the flux tube and its surroundings. As we shall discuss in 
Sect. 6.1, this has important consequences. 

6. Nature of the Evershed Flow 

Previous observational analyses, which neglected the fact 
that in our resolution element we have mixed signals com¬ 
ing from the flux tubes and their magnetic surroundings, 
obtained that the magnetic field strength showed a strong 
radial decrease from Hi nner — 2500 G to H oute r — 1000 
G (Beckers & Schroter 1969; Wittmann 1974; Lites & 
























10 


Borrero et al: On the nature of the Evershed flow 


Skumanich 1990; McPherson et al. 1992; Lites et al. 1993; 
Keppens & Martinez Pillet 1996; Stanchfield et al. 1997). 
This implies that the magnetic field is larger at the inner 
footpoint of a siphon flow carrying loop than in the outer 
one, making it difficult for this mechanism to work. 

If the fine structure of the penumbra is taken into ac¬ 
count (Riiedi et al. 1998, 1999, Borrero et al. 2004) a more 
favourable situation appears: an inclined magnetic field 
whose strength rapidly decreases with radial distance and 
an almost horizontal magnetic field carrying the Evershed 
flow with a far smaller radial drop in field strength are in¬ 
ferred (see Fig. 6 ; bottom panels). The small decrease in 
the field strength for the component carrying the Evershed 
flow still implies a stronger magnetic field at small r/R , 
but now the difference between inner and outer penumbra 
is significantly reduced to ~ 300 — 500 G. Although this 
finding on its own is insufficient to produce an outwards 
accelerated flow it represents a more favourable situation 
for the siphon flow mechanism, specially if we take into ac¬ 
count that, as noticed by Montesinos & Thomas (1997), 
the inner and outer footpoints, as measured by observa¬ 
tions, are not necessarily at exactly the same height. 

6.1. Gas pressure gradient 

According to Fig. 6 (bottom panels) the magnetic field 
strength of the flux tubes is much smaller than that of 
the surrounding atmosphere in the inner penumbra at a 
geometrical height z = Zq (central position of the tube). 
By requiring total pressure balance between the flux tube 
interior and the external atmosphere (Eq. 4) one can eas¬ 
ily deduce that -P gas ,t(2o) ^ Pgas,s(zo)- At large radial 
distances the situation is such that the magnetic field 
strength in the flux tube and its surrounding atmosphere 
becomes very similar and therefore: P gas ,t(zo) — Pgas,s(zo)- 
In Fig. 8 we plot the actual AP gas = P gas ,t(A)) - Pgas, s(~o) 
as a function of radial distance in the penumbra (top panel 
for p = 0.91 and bottom panel for p = 0.51). As can 
be seen that this difference decreases almost linearly with 
r/R , reaching values close to zero near the outer penurn- 
bral boundary. 

Unfortunately, this effect does not imply a radial de¬ 
crease in the gas pressure inside the flux tube, because 
of the unknown radial behaviour of the external pressure 
P g as,s(~o)- In particular, its calculation is ill-posed by the 
boundary condition described in Sect. 2, where a value 
of the gas pressure at the highest photospheric layers is 
prescribed for the surrounding atmosphere. This is done 
for all pixels independently of their position on the spot: 
P g as,s(7ma,x) = Po- Consequently, the geometrical height 
scale z s (t s ) in the surrounding atmosphere does not take 
into account pixel to pixel variations of the absolute ge¬ 
ometrical height scale (i.e. Wilson depression). The ra¬ 
dial variation of the difference in gas pressure between the 
flux tube and the surrounding atmosphere, however, only 
depends on the difference in the magnetic field strength 
between them and is calculated under the condition of 


total pressure balance (magnetohydrostatic constraint). 
Although it cannot directly prove that there is a radial 
decrease in the gas pressure along the flux tube axis, it 
does provide a strong indication that this is indeed the 
case. 

6.2. Cooling flux tubes 

Our results show that in the inner penumbra the flux tubes 
are hotter (by roughly 500-1000 K) than their surround¬ 
ings. As we move to larger radial distances the tempera¬ 
ture becomes very similar inside and outside the flux tube. 
This is in very good agreement with the theoretical predic¬ 
tions of the moving flux tube model (Schlichenmaier et al. 
1998a, 1998b). The result that at a given point, r/Rcx 0.7, 
the flux tube becomes cooler than its surroundings can be 
explained by assuming that the flow becomes adiabatic 
due to a less efficient radiative exchange in the higher parts 
of the anchored flux tube (Montesinos & Thomas 1997). 
The temperature increase seen at larger radial distances 
is explained in terms of a shock front (see Sect. 6.3). 

This picture indicates that at the umbral-penumbral 
boundary, a first set of hot flux tubes (small ones or a 
small number of them according to the tiny filling factor 
inferred; see Fig. 5 bottom panels) , appear. Remarkably, 
these flux tubes that have just emerged already carry 
strong velocities (see Fig. 5 middle panels), perhaps as 
an indication that most part of the plasma acceleration 
has already occurred. Again this is in very good agreement 
with the moving flux tube simulations (see Schlichenmaier 
et al. 1998a, 1998b) where the acceleration takes place in 
the inner penumbra, where hot flux tubes carrying plasma 
at about 8000-10000 K rapidly cool down. The opacity in 
such extremely hot tubes would be too large to make them 
observable; therefore what we see is perhaps, the latest 
stages of this cooling process. Needless to say, hotter flux 
tubes in the inner penumbra were already expected con¬ 
sidering that at small radial distances, the magnetic field 
strength in the flux tubes is much smaller than in the sur¬ 
rounding atmosphere, and therefore, in order to keep the 
horizontal total pressure balance between these two atmo¬ 
spheres the gas pressure in the flux tube component has 
to increase. 

6.3. Shock front 

Another point of special interest is the magnitude of the 
Evershed flow. Theoretical models predict different val¬ 
ues for the speed of the flow inside penumbral flux tubes. 
Stationary siphon flow models (Thomas & Montesinos 
1991,1993; Montesinos & Thomas 1997) distinguish be¬ 
tween subcritical and supercritical velocities with respect 
to the characteristic critical tube speed (Thomas 1988, 
Ferriz Mas 1988). Transitions between both regimes can be 
present within the flux tube, leading to shock fronts. Time- 
dependent simulations of thin flux tubes in the penum¬ 
bra predict supercritical flows in most of the penumbra 
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regardless of whether the flux tubes remain horizontal 
beyond the visible limit of the sunspot (Schlichenmaier 
et al. 1998a, 1998b) or bend back within the penumbra 
(Schlichenmaier 2002). Taking into account the fine struc¬ 
ture of the penumbra, large velocities (> 4 km s -1 ) are 
favoured (Wiehr 1995, del Toro Iniesta et al. 2001; Bellot 
Rubio et al. 2003; Bellot Rubio et al. 2004). Most of these 
values have been obtained under the assumption that the 
magnetic field and the velocity vectors are mutually par¬ 
allel (see Bellot Rubio et al. 2003). 

The absolute flux tube velocity is: Vt = vi os ,t/ cosyt, 
where 7 t is the inclination of the tube’s magnetic field 
vector with respect to the observer (both i>i os ,t and 7 t are 
obtained from the inversion), vt is plotted in Fig. 9 (top 
left panel) for the sunspot at an heliocentric angle of p = 
0.51 4 . Also plotted are the local sound speed c s ~ T 1//2 and 
the tube’s critical speed c t = c s c a / \/c 2 + c a , with c a ~ 
Bp -1 / 2 being the Alfven speed. The velocity in the flux 
tube always remains subsonic, although we cannot rule 
out the possibility for supersonic values in the penumbra 
(del Toro Iniesta et al. 2001; Penn et al. 2003) given the 
limited number of radial cuts we are considering. 

At almost all radial distances, Vt < Ct as well, ex¬ 
cept for a few regions in the inner (r/R cs 0.4), and the 
outer penumbra ( r/R ~ 0.78). Given the error bars, we 
do not deem the first case to be reliable. However, for 
large radial distances (arrow in Fig. 9; top left panel) 
it seems plausible that the velocity becomes supercriti¬ 
cal. Remarkably, after this happens, the velocity suddenly 
drops again by roughly 2 km s^ 1 to subcritical values at 
larger distances: r/R > 0.85. This is accompanied by an 
increase in the flux tube temperature of about 300-400 K 
at these locations (see arrow in Fig. 5; top right panel). 
Indeed, this is the behaviour expected from a transition 
between supercritical to subcritical velocities in penurn- 
bral flux tubes (Montesinos & Thomas 1997) produced by 
a shock front that dissipates kinetic energy and heats the 
gas. The presence of shocks in the siphon flows has been 
predicted by numerous authors (Meyer & Schmidt 1968; 
Degenhardt 1989,1991; Thomas & Montesinos 1991, 1993; 
Montesinos & Thomas 1997; Schlichenmaier et al. 1998a, 
1998b; Schlichenmaier 2002), but this is the first time that 
direct observational evidence of such a shock front within 
the penumbra is provided. 

Shock fronts, if present, are likely to produce an en¬ 
hancement both in the line width (e.g. FWHM) and in 
the equivalent width (see Degenhardt et al. 1993; Solanki 
et al. 1996). In the penumbra, the radial variation of the 
line width is dominated by the Zeeman splitting (magnetic 
field). However, the magnetic held affects the equivalent 
widths of the comparatively weak (i.e. unsaturated) lines 
considered here only very slightly. In Fig. 9 (top right pan¬ 
els) we plot the radial variation of the equivalent width for 
Fe I 15648.5 A. It clearly shows an enhancement (~ 10%) 

4 Results for p = 0.91 are not considered, since dividing by 
the cosine of inclination angles close to 90° enhances any error 
in line-of-sight velocities 


at roughly the same radial distance where the how speed 
becomes supercritical (see vertical arrow). Three different 
radial positions, corresponding to locations before, dur¬ 
ing and after the shock have been marked with open cir¬ 
cles. Intensity prohles for Fe I 15648.5 A at those loca¬ 
tions have been extracted and plotted together in Fig. 9 
(bottom left panel). The intensity profile before the shock 
(dashed line) shows an enhancement in the red wing which 
is produced by the strongly red shifted hux tube contri¬ 
bution (i.e. satellite line; see Stellmacher & Wiehr 1980; 
Wiehr et al. 1986; Wiehr 1995,1997). During the shock 
(solid line) the equivalent width of the redshifted com¬ 
ponent is enhanced. This broadening is likely to be pro¬ 
duced by a new structure in our resolution element that 
our model does not account for. Therefore, the only way 
for the inversion code to make the profile broader (within 
the constraints of the chosen model) is to increase the 
micro and macroturbulent velocities. In Fig. 9 (bottom 
right panel) the line broadening velocity in the flux tube, 

defined as: Abroad = yjv % c t + r^ac.t is plotted. It shows 
a peak exactly at the radial positions at which the flow 
speed becomes supercritical. This provides a strong indi¬ 
cation that this spectral signature really corresponds to a 
shock front in the flux tube at large radial distances in the 
penumbra. 

In our observations the shock seems to occupy 1.2-1.4 
arc sec. The fact that we do not observe a jump over a 
smaller radial range is likely to be caused by smearing ef¬ 
fects introduced by seeing. Finally, it is important to recall 
that this result is based on some assumptions that must 
be considered carefully. In particular the density inside the 
flux tube, pt , which is obtained through the gas pressure 
and temperature by applying the ideal gas equation (see 
Sect. 2), does not satisfy vertical hydrostatic equilibrium. 
The Alfven speed, and therefore the tube’s critical speed, 
are affected by any inaccuracies in the density, which in¬ 
troduces some uncertainties. 

7. Summary and conclusions 

We have presented the first full inversions of infrared 
penumbral spectropolarimetric data with a model that 
accounts for the vertical and horizontal inhomogeneities 
of the penumbral fine structure. This model is a slightly 
modified version of the uncombed model of Solanki & 
Montavon (1993; cf Martinez Pillet 2000) and allows for 
the presence of a randomly orientated flux tube embedded 
in a surrounding atmosphere that harbours a strong mag¬ 
netic field. The main advantages of this model, as com¬ 
pared to those used in Paper I (see also Mathew et al. 
2003; Bellot Rubio et al. 2003) are: a -it contains two dif¬ 
ferent atmospheres in the direction perpendicular to the 
observer (flux tube and magnetic surrounding) which al¬ 
lows the observed multilobed Stokes V profiles to be eas¬ 
ily reproduced; b -these two atmospheres are also present 
in the direction parallel to the observer line-of-sight and 
therefore the model includes gradients along the line of 
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Fig. 9. Top left panel: Radial variation of the flow velocity inside the flux tube (solid line). Local sound speed and 
tube’s critical speed are also plotted (dashed and dashed-dotted lines respectively). The vertical arrow marks the 
position where the flow speed becomes supercritical, r/R — 0.78. Top right panel: radial variation of the equivalent 
width for Fe I 15648.5 A line. The open circles are three selected radial positions before, during and after the shock 
occurs. Bottom left panels: Intensity profiles corresponding to the 3 selected radial positions (top right panel in this 
Figure) before (dashed line), during (solid line) and after the shock (dotted line). Bottom right panels: radial variation 
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sight (in the form of sharp jumps in the physical quantities 
at the flux tube’s boundaries) needed to produce asym¬ 
metric Stokes V profiles {8A A 0); c -flux tubes and their 
magnetic surroundings are coupled to each other using to¬ 
tal pressure balance at all heights; d -the current model 
requires fewer free parameters to reproduce the data. 

Feature b makes this model suitable to investigate the 
vertical size of the penumbral filaments, since the radius of 
the flux tube Rt , which is ultimately linked to the amount 
of area asymmetry generated, is also obtained from the 
inversion. The small amount of 8A shown by the Fe I 
lines at 1.56 /./in is not enough to constrain the position of 
the flux tubes’ upper and lower boundaries, except at the 
outer penumbra where 5A is sufficiently large, allowing 
us to set a rather accurate position for the lower bound¬ 
ary of the flux tube at around logr 5 G [—0.5,0], while 
the upper boundary remains undetected. Although for in¬ 
dividual Stokes V profiles the area asymmetry produced 
by the uncombed model matches the observed one only 
relatively inaccurately, general trends such as the radial 
behaviour of 5A are reproduced fairly well. We point out, 


in agreement with Paper I, that the use of lines that show 
a much larger 5A should help in this matter (see Borrero 
et al., in preparation). 

Feature c has allowed us to detect, for the first time, a 
strong radial decrease in the pressure difference between 
the flux tube and its surroundings that is likely to induce 
an outward directed flow. This is to our mind the strongest 
evidence so far supporting siphon flows (Meyer & Schmidt 
1968; Montesinos & Thomas 1993,1997) as the physical 
mechanism driving the Evershed effect. 

In addition, we have seen that the Evershed flow al¬ 
ready carries velocities as large as v ~ 4 km s -1 in the 
inner penumbra. This has passed unnoticed in previous 
works (see Schmidt & Schlichenmaier 2000; Tritschler et 
al. 2004) where only Stokes I was considered. A possi¬ 
ble explanation is found in the small filling factor and 
inclinations, larger than 90 deg with respect to the ob¬ 
server, of the magnetic field vector in the flow channels at 
the inner penumbra. This geometry produces small line 
asymmetries in Stokes J, but large Stokes V zero crossing 
shifts. At the inner penumbra these fast flows are also as- 
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sociated with hot gas and flux tubes that are somewhat 
inclined with respect to the horizontal, in agreement with 
Schmidt & Schlichenmaier (2000) and Rimmele (2004). 
These results are in close agreement with the dynamical 
simulations of penumbral flux tubes by Schlichenmaier et 
al. (1998a,1998b) and Schlichenmaier (2002). The tubes 
reach the same temperature as their surroundings very 
rapidly in the radial direction, and at the same time the 
flow speed increases smoothly (up to v ~ 5 km s -1 ) as the 
pressure drops (Montesinos & Thomas 1993, 1997). 

At large radial distances the flow speed suffers a sud¬ 
den decrease that is associated with positions where the 
tubes return to the solar interior (Westendorp Plaza et 
al. 1997). In addition to this well established result, we 
have also detected a possible transition between critical 
and subcritical velocities (as predicted by Montesinos & 
Thomas) that is co-spatial with a rise in temperature and 
equivalent width at the outer penumbra, and seems to in¬ 
dicate that part of the kinetic energy is being dissipated 
into thermal energy. Therefore, part of the sudden drop in 
the velocity at the outer penumbra could be ascribed to 
the development of shock fronts. 
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Fig. 7. Top panel: area asymmetry from the best-fit 
Stokes V profiles, SApn, versus area asymmetry of the 
observed Stokes V, SAobSj for both Fe I lines for the case 
/i = 0.91. For /i = 0.51 (not shown) the result is very sim¬ 
ilar. Filled circles are for Fe I 15648.5 A, open circles for 
15652.8 A. Middle panel: average position, in the optical 
depth scale, of the lower boundary of the flux tube (solid 
line). The logri .56 = 0 level is indicated by the dashed 
line. Maximum and minimum individual deviations of the 
10 radial cuts considered are indicated by the shaded area. 
Bottom panel: radial variation of the fitted (dashed line) 
and observed area asymmetry (solid line) of Fe I 15648.5 

A. 
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fi= 0.91 




Fig. 8. Radial variation of the gas pressure difference, 
AP gas , between the flux tube and its surrounding atmo¬ 
sphere at 2 = zq. Top panel: for NOAA 8706 at an he¬ 
liocentric angle of /.t = 0.91. Bottom panel: the same for 
/t =0.51. 














